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Abstract 

We  study  the  acoustic  scattering  properties  of  a  phononic  crystal  designed  to  behave  as  a  gradient 
index  lens  in  water,  both  experimentally  and  theoretically.  The  gradient  index  lens  is  designed 
using  a  square  lattice  of  stainless-steel  cylinders  based  on  a  multiple  scattering  approach  in  the 
homogenization  limit.  We  experimentally  demonstrate  that  the  lens  follows  the  graded  index 
equations  derived  for  optics  by  mapping  the  pressure  intensity  generated  from  a  spherical  source  at 
20  kHz.  We  find  good  agreement  between  the  experimental  result  and  theoretical  modeling  based 
on  multiple  scattering  theory. 

PACS  numbers:  43.20.Fn,  43.58.Ls,43.20.Dk 
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Composed  of  ordered  arrays  of  scatterers  similar  to  atoms  in  a  conventional  solid, 
phononic  crystals  (PnC)  are  a  class  of  metamaterial  designed  to  control  acoustic  wave 
propagation  in  a  medium.  PnCs  have  been  proposed  for  a  broad  range  of  applications 
in  wave  acoustics,  with  acoustic  lensing  [M2]  featuring  prominently  in  the  literature  due 
in  part  to  the  ease  with  which  focusing  can  be  achieved  by  altering  a  crystal’s  gemoetric 
shape  BE!  or  compositional  structure.  HH6]  Although  negative  index  lenses  have  received 
much  attention  due  to  their  potential  for  near-held  imaging,  HH9]  some  positive  index  solu¬ 
tions  such  as  the  acoustic  analogue  of  the  optical  graded  index  lens  mmm  have  not  yet 
been  explored  experimentally.  In  addition,  the  majority  of  PnC  experiments  have  been  per¬ 
formed  in  air,  PO  El  H21 H3]  where  the  large  density  contrast  with  respect  to  the  constituent 
scatterers  in  the  crystal  (typically  metals)  allows  the  scatterers  to  be  treated  as  rigid.  We 
demonstrate  below  that  despite  the  physical  limitation  in  impedance  contrasts  between  an 
aqueous  medium  and  the  scattering  elements,  it  is  possible  to  successfully  design  a  PnC 
that  behaves  as  an  ideal  graded  index  lens  (GIL)  in  water  based  on  a  fully  elastic  multiple 
scattering  theory  (MST).  [U  [CL  ;T4] 

Figure  [jja)  shows  a  plan  view  schematic  of  the  GIL  design.  The  axes  of  Fig.  [l](a)  and 
throughout  the  paper  are  oriented  with  the  lens  center  at  position  (x,y)  =  (0,0).  The  GIL 
is  made  up  of  75  stainless  steel  cylinders  (T-316)  that  are  75  cm  in  length  and  arranged  in 
a  square  lattice  with  spacing  a  =  1.8  cm  and  dimensions  5a  x  15a.  Figure  [TJ(b)  plots  the 
cylinder  radii  R{y),  which  are  stepped  toward  zero  at  each  successive  layer  above  and  below 
the  central  axis  ( y  =  0)  of  the  GIL.  In  the  homogenization  limit  (propagation  wavelength 
A  >  4a)  each  stratified  layer  in  Fig.  [jja)  can  be  treated  as  an  effective  medium.  MST  [4] 
is  used  to  calculate  each  layer’s  effective  sound  speed  ce//,  which  is  inversely  proportional 
to  the  filling  fraction  of  the  cylinders.  The  layers  will  have  an  effective  refractive  index 
neff  =  Cb/ceff  ( Cb  =  1470  m/s  is  the  sound  speed  in  water)  that  is  maximal  at  the  center  of 
the  GIL  and  decreases  to  that  of  water  at  the  edges.  Our  choice  of  R{y)  in  Fig.  [ljb)  produces 
a  graded  ne//  that  obeys  the  same  relation  as  an  optical  GIL,  [15]  ne//  =  no(l  —  a2y2)1'/2, 
where  no  is  the  refractive  index  at  the  central  layer.  Our  design  results  in  no  =  1.2  and 
a  =  0.04  cm-1. 

Figure  [l](c)  shows  an  image  of  the  GIL  submerged  in  a  6  x  6  x  4  m3  isolation  tank.  The 
cylinders  are  mounted  between  reinforced  Plexiglas  plates  to  provide  stability.  Acoustic 
waves  are  produced  by  a  10  cm-diameter  spherical  source  at  20  kHz  (A  ~  4a).  The  wave 
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FIG.  1.  (a)  Plan  schematic  of  the  gradient  index  lens,  (b)  Cylinder  radius  R  plotted  vs  position 
along  the  y-axis.  (c)  Digital  photograph  of  the  GIL  in  the  isolation  tank. 


propagation  is  measured  in  the  time-domain  using  monitoring  hydrophones  at  a  sampling 
rate  of  1  MHz.  Hydrophones  are  mounted  to  the  source  and  onto  a  translational  3-axis 
Velmex  VXM®  positioning  system.  The  transmission  intensity  is  measured  by  averaging 
over  a  10-cycle  pulse  from  the  source;  this  pulse  is  long  enough  to  approximate  a  continuous 
wave  measurement,  while  being  short  enough  to  prevent  contamination  from  reflections  off 
the  surfaces  of  the  tank.  We  have  experimentally  verified  that  the  intensity  P0  produced  by 
the  source  in  the  absence  of  the  GIL  drops  radially  in  proportion  to  1/r2. 

Figure  [2|a)  shows  the  normalized  pressure  amplitude  P/Pq  measured  after  transmission 
through  the  GIL  on  the  side  opposite  the  source  {x  >  0).  The  source  is  located  at  (x,y)  = 
(—196,0)  cm,  and  both  the  source  and  the  translational  hydrophone  are  positioned  in  the 
plane  bisecting  the  axial  center  of  the  cylinders.  The  GIL  is  shown  schematically  to  scale  and 
at  its  proper  location  in  each  figure  throughout  the  paper.  The  data  in  Fig.  |2](a)  is  measured 
2.144  ms  after  the  first  cycle  began  to  leave  the  source.  This  time  gives  a  snapshot  when 
the  pulse  is  centered  on  an  enhancement  in  signal  amplitude  observed  in  the  vicinity  of 
x  ~  80  cm. 

Figure  [2](b)  shows  the  normalized  intensity  averaged  over  the  10-cycle  pulse  and  obtained 
from  the  same  data  set  shown  in  Fig.  |2](a).  As  in  Fig.  |2](a),  a  clear  focusing  peak  is  observed 
centered  close  to  x  ~  80  cm.  In  Fig.  [2|c)  we  show  a  two-dimensional  MST  calculation  ill 
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FIG.  2.  (a)  Normalized  pressure  amplitude  P/ Pq  plotted  vs  x  and  y ,  measured  2.144  ms  after  the 
initial  pulse  leaves  the  source,  (b)  Measured,  normalized  pressure  intensity  ( P/Pq )2  plotted  vs  x 
and  y  after  averaging  over  a  10-cycle  pulse,  (c)  Normalized  pressure  intensity  ( P/Pq )2  calculated 
using  multiple  scattering  theory.  The  focusing  peak  maximum  is  marked  by  a  +. 


of  the  total  pressure  intensity  (incident  plus  scattered)  derived  by  placing  a  continuous- wave 
cylindrical  source  at  (x,y)  =  (—196,0)  cm.  The  calculation  assumes  the  cylinders  to  be  a 
penetrable  elastic.  mm  As  with  the  experimental  data,  the  simulated  pressure  intensity  is 
normalized  to  that  of  the  source  in  the  absence  of  the  GIL.  The  source  amplitude  is  a  Hankel 
function  P0  =  H^\kr)  with  wavevector  k  =  n/2a.  The  MST  simulation  also  shows  a  clear 
focusing  peak,  but  with  two  important  differences:  (1)  the  measured  intensity  is  ~  2  times 
larger  than  that  of  the  simulation,  and  (2)  the  simulated  focusing  peak  is  slightly  farther 
from  the  GIL  and  decays  more  slowly. 

To  quantify  whether  our  GIL  design  behaves  as  an  ideal  lens,  in  Fig.  [3|a)  we  present 
measurements  of  the  focusing  peak  along  the  central  axis  of  the  lens  ( y  =  0)  for  different 
source  positions  ds.  For  each  ds,  a  large- amplitude  focusing  peak  is  observed  above  x  > 
60  cm,  while  smaller  peaks  are  also  observed  closer  to  the  GIL.  As  the  source  is  moved 
closer  to  the  GIL,  the  large-amplitude  peak  moves  away  in  qualitative  agreement  with  the 
expected  behavior  of  a  lens.  In  Fig.  [3](b)  we  show  MST  calculations  along  y  =  0  for  source 
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FIG.  3.  (a)  Measured,  normalized  pressure  intensity  vs  x  for  source  positions  ds  =  196,  185.8, 

175.7,  165.5,  155.4,  145.2,  135,  124.9,  and  114.7  cm.  (b)  Normalized  pressure  intensity  calculated 
using  multiple  scattering  theory  for  source  positions  ds  =  109a,  103a,  98a,  92a,  86a,  81a,  and 
75 a.  Inset:  expanded  region  of  the  y-axis  showing  the  focusing  peak  positions,  (c)  Two  focusing 
peaks  from  panel  (a)  are  replotted  as  circles  (upper  region  offset  for  clarity).  Black  lines  indicate 
a  double-gaussian  fit,  with  blue  (Peak  1)  and  red  (Peak  2)  lines  showing  the  component  gaussians 
individually,  (d)  Inverse  peak  positions  l/dpi,2  plotted  vs  inverse  source  positions  l/ds.  Blue  and 
red  lines  are  fits  to  the  trends  of  Peaks  1  and  2  respectively.  The  dashed  line  plots  the  peak 
positions  of  the  simulated  data  in  panel  (b). 

positions  similar  to  those  in  the  experiment.  On  first  inspection  it  appears  that  the  theory 
shows  slowly  decaying  focusing  peaks  that  change  very  little  with  ds.  However,  expansion  of 
the  region  around  the  focusing  peaks  [Fig.  |3](b)  inset]  reveals  that  the  peak  positions  move 
away  in  a  manner  similar  to  the  experiment. 

We  now  analyze  the  experimental  data  in  Fig.  |3](a)  above  x  >  62  cm  to  determine  whether 
the  focusing  positions  in  this  region  follow  the  ideal  lens  equations.  For  an  ideal  lens,  the 
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focusing  peak  positions  dp  should  scale  with  ds  as  1  /  dp  =  1/ f  —  l/ds.  The  focal  length  /  of 
a  GIL  can  be  approximated  as,  [TBj 

/« - - -  (1) 

n^a  sm  at 

where  t  =  5a  is  the  thickness  of  the  GIL.  Equation  (JTJ  gives  an  estimate  of  /  =  58.9  cm 
using  the  values  of  n0  and  a  calculated  in  the  effective  medium  approximation.  [TU] 

A  close  inspection  of  Fig.  |3](a)  reveals  that  the  data  above  x  >  62  cm  is  actually  composed 
of  two  superimposed  peaks  that  both  move  to  larger  x  as  a  function  of  ds.  [L7j  Figure  |3](c) 
shows  two  examples  of  a  double-gaussian  fit  to  the  data  in  this  region  using  a  standard 
unconstrained,  nonlinear  optimization  routine.  The  gaussians  resulting  from  the  fit  are 
shown  individually  (blue  and  red)  in  addition  to  the  combined  fit  (black).  Although  the  fit 
equation  71  e~''3'  (x~dr'  f  +7 2e~^x~dP2'>2  contains  six  free  parameters,  the  purpose  of  the  fit  is 
to  obtain  an  estimate  of  the  peak  positions  dp  17  and  their  relative  amplitudes  772.  Fig.  |3](c) 
demonstrates  that  the  data  is  well  described  by  the  double-gaussian,  with  a  low-amplitude 
peak  (Peak  1)  closer  to  the  GIL  and  a  larger-amplitude  peak  (Peak  2)  farther  away.  In 
both  cases  the  amplitude  of  Peak  2  is  about  three  times  larger  than  Peak  1,  suggesting  that 
Peak  2  is  the  main  focusing  peak  of  the  GIL.  The  relative  amplitudes  of  Peaks  1  and  2  are 
observed  to  follow  the  same  qualitative  behavior  for  all  the  source  positions  in  Fig.  [3|a). 

Figure[3](d)  plots  the  inverse  positions  l/dpii2  extracted  from  the  gaussian  fits  as  a  function 
of  l/ds.  An  ideal  lens  will  produce  a  linear  trend  with  a  slope  of  —1  and  an  intercept  of 
1  //.  Although  the  trends  for  both  peaks  are  linear  and  have  intercepts  that  yield  similar 
focal  lengths,  the  slope  of  Peak  1  is  much  less  than  that  of  an  ideal  lens.  However,  the 
trend  for  Peak  2  results  in  a  slope  of  —1  and  its  focal  length  /  =  59.3  ±  1.5  cm  agrees  with 
the  estimate  of  /  calculated  using  Eqn.  ([!]).  The  dashed  line  in  Fig.  |3](d)  plots  the  peak 
locations  obtained  from  the  MST  calculations  in  Fig.  [3](b).  The  theory  produces  a  slope  of 
— 1  and  focal  length  /  =  61.1  cm  that  closely  match  both  the  measured  data  and  the  ideal 
lens  equations. 

We  propose  that  Peak  1  and  the  other  low-amplitude  peaks  in  Fig.  |3](a)  are  the  result 
of  constructive  interference  between  waves  scattered  from  the  support  structure  of  the  lens. 
Low-amplitude,  circular  interference  fringes  can  be  observed  in  Fig.  [2](a)  emanating  from 
above  and  below  the  plot  area  centered  at  x  ~  25  cm.  These  fringes  are  the  result  of 
scattering  off  of  stabilizing  pillars  at  the  corners  of  the  GIL  support  structure.  While  this 
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FIG.  4.  (a)  Measured,  normalized  pressure  intensity  ( P/Pq )2  plotted  vs  x  and  y  with  the  source 

at  a  14.7°  angle  with  respect  to  the  x-axis.  (b)  Normalized  pressure  intensity  ( P/Pq )2  calculated 
using  multiple  scattering  theory  with  the  source  at  a  15  °  angle.  White  lines  indicate  the  off-axis 
angle;  positions  of  the  expected  focusing  peaks  are  marked  with  a  +. 

interference  can  be  reduced  by  limiting  the  averaging  to  only  the  first  few  pulse  cycles,  a  small 
number  of  cycles  gives  a  poor  approximation  to  a  continuous  wave  measurement  and  limits 
the  number  of  multiple  scattering  events  that  contribute  to  the  focusing  peak.  Therefore  we 
have  chosen  to  average  over  many  cycles  and  rely  on  the  gaussian  fitting  routine  to  remove 
the  spurious  interference. 

Figure  [4]  demonstrates  that  our  GIL  design  acts  as  a  lens  with  the  source  off  the  central 
axis.  Figure  |4|a)  shows  (P/Po)2  measured  with  the  spherical  source  located  at  a  14.7°  angle 
with  respect  to  the  origin.  Figure  [4^b)  shows  the  MST  calculation  for  the  same  source 
location.  Thin  white  lines  in  Figs.  |4](a,b)  indicate  a  14.7°  angle  with  respect  to  the  x-axis 
and  extend  to  the  expected  focusing  positions  of  an  ideal  lens  with  /  «  60  cm.  Both 
the  measured  data  and  the  MST  calculation  demonstrate  a  strong  focusing  peak  at  the 
expected  location.  Note  that  interference  fringes  from  the  support  pillar  are  clearly  evident 
superimposed  on  the  focusing  peak  in  Fig.  |4|a). 

In  summary,  we  have  designed  and  constructed  a  gradient  index  lens  that  operates  in 
water  at  sonic  frequencies.  Our  transmission  measurements  demonstrate  that  our  GIL  design 
focuses  as  an  ideal  lens  based  on  the  optical  GIL  equations.  Our  measurements  are  also 
consistent  with  the  focusing  positions  obtained  from  two-dimensional  models  using  multiple 
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scattering  theory.  We  emphasize  that  our  GIL  behaves  as  an  ideal  lens  at  the  limit  of 
homogenization  (A  ~  4a)  and  with  a  thickness  on  the  order  of  a  wavelength  (f  =  5A/4). 
Such  performance  at  the  limit  of  homogenization  theory  demonstrates  the  versatility  of 
phononic  crystals  designed  using  multiple  scattering  theory. 

This  work  was  supported  by  the  U.S.  Office  of  Naval  Research. 
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